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in"' Abstract 

■ We study resummation of transverse-momentum-related large logarithms generated from soft- 

. gluon radiations in soft-collinear effective field theory. The anomalous dimensions of the effective 

quark and gluon currents, an important ingredient for the resummation, are calculated to two-loop 
order. The result at next-to-leading-log reproduces that obtained using the standard method for 
ly^ ' deep-inelastic scattering, Drell-Yan process, and Higgs production through gluon-gluon fusion. We 

comment on the extension of the calculation to next-to-next-to-leading logarithms. 
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I. INTRODUCTION 



Perturbative calculations of high-energy processes with widely separated scales often yield 
large logarithms involving the ratios of the scales, which shall be resummed to all orders 
to achieve reliable predictions. In processes like deep- inelastic scattering (DIS), Drell-Yan 
(DY) and the Higgs-boson production, multiple hard scales appear in a certain kinematical 
limit such as threshold production region 0,1^, and/or when a low transverse momentum 
0, 0] of the final states is measured. For the transverse momentum distribution, the rigor- 
ous theoretical study in QCD started with the classical work on semi-inclusive processes in 
e^e~ annihilation by Collins and Soper in where a factorization was proved based on the 
transverse momentum- dependent (TMD) parton distributions and fragmentation functions 
[^. The resummation of TMD large logarithms was performed by solving the relevant en- 
ergy evolution equation. This approach was later applied to the DY process in Q, where a 
general and systematic analysis of the factorization and resummation were performed. This 
latter procedure became known as "CoUins-Soper-Sterman" (CSS) resummation formalism. 
The resummed formulas are used for many processes with the relevant coefficients extracted 
by comparing between the expansion of the resummed expressions and the fixed-order cal- 
culations 0. Although the factorization approach is sound and rigorous, it involves often 
quantities that are not manifestly gauge invariant. Moreover, the physics of the energy 
evolution does not seem transparent. 

In this paper we pursue the same resummation from a different path. We exploit the 
fact that there are (at least) two well- separated hard scales which are naturally appropriate 
for an effective field-theoretic approach to perform the resummation. The recently proposed 
"soft-coUinear-effective-theory" (SCET) |^ is useful here. Although it was originally applied 
for the study of heavy B meson decays, it was later generalized to other high-energy processes 
[lol | . Recently, this effective theory has been used to study the threshold resummation for the 
DIS structure function as x — 1 in and for DY ^^l- The resummation in the effective 
theory is performed by studying the anomalous dimensions 71 of the effective operators after 
performing matching between the full and effective theories. The exponentiated Sudakov 
form factor appears by running down the scale of the matching coefficient from the higher 
scale jji'^j ~ down to the lower scale /i| ~ (A ^ 1). 

A SCET study of transverse-momentum dependence for Higgs-boson production was 
initiated in Refs. where the authors work directly in momentum space, and the next- 
to-leading logarithmic result was deduced from a full QCD calculation. Here we consider the 
same resummation for the Drell-Yan process as well as the standard-model Higgs production 
at hadron colliders. Instead of working in the momentum space, we use the impact parameter 
space. Moreover we perform calculations directly in SCET. 

In general, the resummation formula for these processes can be written in the following 
form, taking the DY process as an example 0, 



da 

2 = ^0 



dQ'^dydQ. 



J (^'''"^"' + 



where ao = Att'^oI^/ (9s(5^) represents the Born-level cross section for DY production. The 
variables used here are standard: is the invariant mass of the DY pair; Qt is the observed 
transverse momentum relative to the beam axis; y = (1/2) ln{Q^ + Q^)/{Q^ — Q^) is the 
rapidity; s = {Pi -1-^2)^ is the center-of-mass energy squared with Pj the momentum of 
the incoming hadrons. The variable Xi and X2 are the equivalent parton fractions, Xi = 
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^ Q'^l s and X2 — y/ Q'^ /s e~^. The W term contains the most singular contributions 
at small Qt, resummed to all orders in perturbation theory. The second term Y{Qt, Q^) 
represents the regular part of a fixed-order calculation for the cross section, which becomes 
important when the transverse momentum Qt is on the order of a/Q^. 

The main result of our study can be summarized by the following formula for the W 
term, 

qq' 

X {Cg ® q) (xi, 6, g, nl) X {Cg> ® q') (x2, 6, g, ijI) , (2) 

where q and q' are the parton distributions and/or fragmentation functions related to the 
processes studied, and the notation stands for convolution. Two matching coefficients 
appear in the above formula: one is C(g^, as{Q'^)) connecting between the matrix element of 
full QCD current and the effective theory analogue at the scale Q^; the other is the coefficient 
function Cg obtained by calculating the processes in SCET at a lower scale /i| ~ Q^- The 
exponential suppression form factor arises from the anomalous dimension of the effective 
current in SCET, 

5(g^/.i)= ^^27l(/.^a.(/.^)) (3) 

The anomalous dimension is identical for DIS and DY processes, and depends only on the 
effective theory operators. Moreover, the same 71 controls the threshold and the low trans- 
verse momentum resummations. The matching coefficient C(g^, as{Q'^)), on the other hand, 
is process dependent (but independent of threshold or transverse-momentum resummation) . 
All the large double logarithms are included in the Sudakov form factor S. 

In the main body of the paper, we show how to get the above resummation formula using 
SCET. In Sec. 11, we consider the resummation for the DY process in SCET. The anomalous 
dimensions of the effective current are calculated up to two-loop order. A comparison will 
be made with the CSS formahsm. In Sec. Ill we will briefly discuss the extension of the 
formalism to the standard-model Higgs production. We conclude the paper in Sec. IV. 

II. DRELL-YAN PRODUCTION AT LOW Qt IN SCET 

In soft-collinear effective field theory, we consider Drell-Yan production with finite trans- 
verse momentum Qt in two steps, assuming Q 3> Qt 3> Aqcd- In the first step, one 
integrates out all loops with virtuahty of order Q to get an effective theory called SCETj in 
which there are only coUinear and soft modes. The coUinear modes have virtuality of order 
Qt- In the second step, one integrates out the coUinear modes with virtuality Qt- In this 
case, the theory is matched onto SCETn which is just the ordinary QCD without external 
hard scales. The soft physics is now included in the parton distribution functions. 

Let us consider the first step: integrating out modes of virtuality of order Q. At low 
Qt, the most singular contribution in DY comes from the form-factor type of diagrams, in 
which the quark and antiquark first radiate soft gluons, followed by an annihilation vertex 
decorated with loop corrections. If the gluon radiations are attached to the loops, the soft 
gluon limit does not give rise to any infrared singularity, and hence diagrams yield higher- 
order contributions in Qt/Q- Therefore, one needs to consider only the form-factor type 
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of diagrams in studying the virtual corrections. To integrate out the hard modes (where 
all gluon momenta are of order Q- see, e.g., 0] ) from the theory, we can match the full 
QCD current onto the (gauge invariant) SCET current 3| and the matching coefficient 
contains the hard contributions. By exploiting the non-renormalizability of the form factor, 
we write down a simple renormalization group equation from which we can extract the 
anomalous dimension of the effective current. 

In the second step, we calculate the SCET cross section. One needs to compute only 
the real contributions since the virtual diagrams are scaleless in the effective theory and 
vanish in on-shell dimensional regularization (DR). At this stage the cross section obtained 
is the same as the one in full QCD taken to the relevant kinematical limit. This has been 
established at one-loop order in the threshold resummation for DIS and DY jlll . This 
will also be the case for DY in low transverse momentum limit. 

Resumming the large logarithmic ratios is performed by considering the scale depen- 
dence of the matching coefficient, which is controlled by the anomalous dimension of the 
effective current. By running down the scale from ~ to the lower scale all 
the large logarithms exponentiate and give rise to the Sudakov form factor. In the follow- 
ing, we will demonstrate how this can be done systematically, providing a powerful tool 
for future studies. The basic lagrangian and Feynman rules for SCET can be found in 
Refs.0. In our calculation, we choose the two light-like vectors: n = l/\/2(l,0,0, —1) and 
n = l/\/2(l, 0, 0, 1). The matching is made between the full QCD current ip'j^ip and the 
SCETi current ^nWn'j^W^C.n, where C,n,n are the collinear quark fields in SCET, and Wn^ 
are the collinear Wilson lines. The collinear Wilson lines appear as a requirement of the 
collinear gauge invariance in the SCETj lagrangian. All calculations are performed in Feyn- 
man gauge and MS scheme in c? = 4 — 2e. DR is employed to regularize both infrared (IR) 
as well as ultraviolet (UV) divergences. 



A. Matching at scale and anomalous dimension of SCET current at two loops 



At the scale of Q^, we need to consider only the virtual contributions to the quark form 
factor both in full QCD and the effective theory. The full QCD calculation of the quark form 
factor has so far been done up to two-loop order . On the other hand, the virtual diagrams 
in the effective theory are scaleless, and the relevant form factors vanish in DR where the 
UV and IR divergences cancel out at every order in a^. Since the effective theory captures 
the IR behavior of the full theory (and the UV divergences are cancelled in both theories by 
respective counterterms) , the matching coefficients for the effective theory operators will be 
obtained from the finite parts of the full QCD calculation. In the following, we will study 
the matching coefficients at one- and two-loop orders. 

At one- loop order, the matching coefficient has been calculated (C = X^n^*""^) IH) 3, 



2 31n — 



6 



(4) 



where, throughout the paper, is the renormalized running coupling constant. 

From the above result, we can get the order-a^ anomalous dimension for the SCET current 



(Ti = E„7!"') 



7!" 



TT 



^ 2 



(5) 
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The matching coefficient satisfies the renormahzation group equation, 



dln/i 



At the scale /x^ = we have, 

C(Q^a,(Q2)) = l + C(^) = l 



2tt 
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(6) 



(7) 



We note that the matching scale can also be chosen as /xf^ = ClQ^ with C2 a constant of order 
unity. The C2 dependence of the matching coefficient corresponds to the C2 dependence in 
the CSS resummation 0]. However, in order to minimize the logarithms in the matching 
coefficient C(C|(5^, as(C|Q^)), the best choice seems to be C2 = 1. 

To get the anomalous dimension at two-loop order, we need to calculate the matching 
coefficient up to the logarithmic term at the same order. From the quark form factor at 
two-loop order, we get the matching coefficient for DY process. 



(4vr) 



_!|__31n-^^ + 8-C(2) 



+ ( 24C(3) + ^ - 12C(2) 



+NfCp 

+CaCf 

_ / 254:5 
~ [~54' 



In 



In^ 



2-Q 



11, 3 

— In"* 



9 



+ 2 ( 4C(2) - 



209 4 



332 



27 3 



C(2) In 



In^ 



+ 



22 



C(2) - 26C(3)) In^ 



(8) 



Here ■ ■ ■ stands for the constant terms which do not contribute to the anomalous dimen- 
sion. The result contains an imaginary part because the DY form factor contains final-state 
rescattering. From the above result we can calculate the anomalous dimension at order a^. 
The expansion of Eq. (6) to order gives 



9 .^(2)^ 9a, d 



d\nn dlnfida 
The running of the strong coupling constant Q;s(/i) is 

d 



CW ^ ^(1)^(1) + ^(2)c'{0)= ^(1)^(1) + 



7i 



7i 



7i 



(2) 



where [3, has the following expansion. 



a. 



2 

S 1 



Air 



a. 



(3q and (3i are defined as 
11 



/3o 



34 



20 

Cl - ANfTpCp - —NfTpCA 



(9) 



(10) 



(11) 



(12) 
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where Nj is the number of hght flavors, Ca = Nc with Nc the number of colors, Cp = 
{Nq — l)/2Nc, Tp = 1/2. Using the above formulas, we get the anomalous dimension at 
two-loop order. 



(2) 

7i 



67 

36" I2J 



-Nf 



/i2 



+ 



13 



C(3) 



961 



11 



16 X 27 4^ 



TT^ CaC^ 



TT 

+ I 777 + 



65 



+ 



24 8 X 27 

TT 



NfC 



^2 3 



4 16 

Using the notation of [16|, we can write 



3C(3) 



7i 



(2) 



a: 



1 (1) , (1) 



(47r)2 V2 '^5" "2 



2/, 



g2 



where 



IKqq 



16CpK, K = Ca(^^- C(2)) - ^NfTp 



is proportional to the anomalous dimension of a Wilson line cusp, and 

17 8 



.(1) 



Tqq' = (3-24C(2)+48C(3)) + C^C^ 
+NfTpCF (-^ - |C(2)) 



C(2) - 24C(3) 



(13) 



(14) 



(15) 



(16) 



is proportional to the coefficient of 5(1 — x) in the quark splitting function, and finally 



(1) 

q2 



404 

-2/3oC(2) + Ca[ — - 14C(3) ) + NjTpCp 



112 
"27 



(17) 



The above calculation can be repeated for the DIS process and the same anomalous 
dimension is obtained. 



B. Matching at scale ^\ 

To perform matching at the lower scale <^ we calculate the cross section in SCETj 
and match it to a product of quark distributions. From the result we can extract the 
coefficient functions. The parton distributions in SCETn are the same as those in full QCD 
11 . 1^ 1 . Below the scale Q^, the hard modes have been integrated out, and have been taken 



into account by the matching condition at . Therefore, the calculation of the cross section 
at /x| is performed with SCETi diagrams, including both virtual and real contributions. As 
mentioned earlier the virtual diagrams in SCET are scaleless and vanish in pure DR. As 
such one can ignore them, but the counterterm for the effective current must be taken 
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(a) (b) (c) id) (e) 

FIG. 1: Non-vanishing Feynman diagrams contributing to Drell-Yan production in the soft- 
coUinear- effective theory: (a) for the soft gluon radiation; (h)-(e) for n and n coUinear gluon 
radiations. The mirror diagrams of (a-c) are not shown here hut are included in the results. 



into account. This is equivalent to taking the UV-subtracted contribution from the virtual 



diagrams. The counterterm has been calculated in jll 



c.t. 



1 3 1 /"^ 



2e 



(18) 



where /x^ will be taken as 

The real contribution contains coUinear and soft gluon radiation diagrams. Some of these 
diagrams are identically zero because of = or = 0, or because of the equation of 
motion for the external coUinear (anti-)quark. The non-vanishing diagrams are shown in 
Fig. 1, including soft-gluon interference contribution, n and n coUinear gluon radiations 
and their interferences. The contribution from Fig. 1(a) is, including that from the mirror 
diagram. 



(fQxdy 



fig.l(a) — CTQ" 



agCp 1 



Ql 



Xi] 



(19) 



Fig. 1(b) represents the interference between the n-coUinear gluon radiation with coUinear 
expansion of the current operator, and its contribution is. 



|fig.i{b) — ctq- 



Qi 



Xi6{x2 "~ 1) _|_ ^{^1 ~ 1) 



Xi 



1-X2) 



6{xi - 1)6{X2 - 1) In 



(20) 



while Fig. 1(c) for the ri-coUinear gluon radiation yields. 



d^Qxdy 



|fig.l{c) — O'O" 



Qi 



X2S{Xi — 1) 6{X2 — 1) 



X2 



l-Xi, 



6{xi - 1)6{X2 - l)ln 



(21) 



Figs. 1(d) and (e) stand for the n and n coUinear gluon radiations, respectively, and their 
sum is, 



d^Qrdy 



I fig. 1 (do) 



(y. C X 

Q"o 9 7^ (1 - e) [(1 - xi)5{x2 - 1) + (1 - X2)5{xi 



27r2 Q 



1)] 



(22) 



T 



The sum of the above contributions reproduces the result of the full QCD calculation in the 
limit of low transverse momentum (see, e.g., [3l)- 

From the above, the real contribution contains soft divergences (i.e, when Q\ 0), which 
will be cancelled by the virtual contribution. To see this cancellation explicitly, we Fourier- 
transform the cross section from the transverse momentum space into the impact parameter 
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6-space. The result is W{b, Q^, in SCET, including real and virtual contributions (the 
contourterm, Eq. (fTHj) ). 



W{b,Q'^,ljl,Xi,X2) 
= 5{xi - l)S{x2 - 1) + 



27r 

asCp 

27T 



[(1 - Xi)6{x2 

5ixi - l)5ix2 - 1 



1) + (1 



[Vq/q{xi)5{x2 - 1) + (xi ^ X2)] (^-^ - 7e + In ^^^25^ ) 



X2)5(xi - 1)] 



In^ 



+ 21n^ln 



_ 9 TT^' 

4 ^ y 



where Pg/g is the one-loop quark splitting function, 



V, 



x] 



X 



X 



(23) 



(24) 



The soft divergences in 14^(6, Q^,//^) have been cancelled. There is, however, the collinear 
divergence left which can be absorbed into the quark distribution at one-loop order. The 
cross section depends on the ultraviolet scale It is somewhat surprising that the above 
result also depends on InQ^, but this is expected from the kinematical constraints of the 
process. 

In order to eliminate the large logarithms in the coefficient function, the best choice for 
the scale fiL is = Cl/b^ with Ci = 2e~'^^. In addition, W{b,Q'^, fijj no longer depends 
on at this order. This, however, may not be true at higher orders because one cannot 
eliminate all InQ^ by making a single choice of Considering the quark distribution at 
one-loop order, we can cast Eq. (j^^ into the form, 

W{b, hI = C^/b\ xi, X2) = Cg{b, xi, fil) ® q{x, fil)Cg{b, X2, fil) O q{x, fil) , (25) 



where Cq reads 



C,{b,x,fil) = 6{l-x) + 



27r 



'l-x)-5{x-l] 



12 



(26) 



and Cq = Cq. Note that we have set the quark charge e„ to 1. 



C. Resummation and comparison with conventional approach 



The final result for W{b, Q^) is a combination of the factors we have calculated in the 
previous subsections, and as have been advertised in the introduction [l^ , 

W{b, Q\ xi, X2) = J2 «.(g'))e-'5(«'' {Cq ® q) (xi, b, Q\ /xi) (Q ® q) (x2, 6, Q^ 



(27) 

where the exponential suppression factor can be calculated from the anomalous dimension 
calculated in Sec. IIA, 



S 



ML 





r 0^ 1 

A\n^ + B 


Li 





(28) 
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Here A and B are the single logarithmic and constant terms in the anomalous dimension 
7i, respectively. The separation of the anomalous dimension into a sum of these two terms 
holds to any order in perturbation theory. The Sudakov form factor comes from the running 
of the matching coefficient C{Q'^,fj,'^) from //^ = down to the scale of /i| which can be 
taken as Cf/b'^. 

The above is the final resummation result in the effective theory. A and B coefficient 

functions can be expanded as a series of a^: A = {^Y and B = X]i=i-^*'*'' {^Y- 

From the result in Sec. IIA, we have the first two terms of these expansions, 



^(1) 
5(1) 



Cf 
3 



5(2) 



67 
36 

?C(3) 



TT 

12 



961 



16 X 27 48 



—n'^ 1 CaC^ 



+ 



71 

T 



16 



3C(3) C| 



71 

24 



65 



8 X 27 



NfC 



(29) 



In the present formulation, the summation of leading logarithms (LL) involves A^^^^; that of 
next-to-leading (NLL) logarithms involves A^"^^ and 5^^); and that of next-to-next-to-leading 
logarithms (NNLL) involves A'^'^\ B^'^\ one-loop Cq, and part of the two-loop Cq [we write 

c, = j:nCi''\ 

To compare with the CSS approach we follow the procedure outlined in [l8| by absorbing 
the factor C{Q'^,as{Q^)) into B and C functions, for example, up to order a^, 



^g/CSs(^) 



Cj^^)(x)+5{l-x)^{C{Q\a.{Q'))Y 



B. 



(2) 



CSS - B^''>+(3o^{C{Q\as{Q'))Y 



(30) 



At two-loop level and beyond, one has to shuffie the In Q^- dependent part of Cq into -Bess 
as well. In this way, we get the CSS resummation formula as. 



e "^^"^ (Qcss ® q) (xi, 6, nl) {Cq/css ® q) {^2, 6, f^l) 



with 



sud 



A 1 

Acss In 



B. 



CSS 



^css will be the same as our A functions, as does B^^\ For C^/css, "we have 



(31) 
(32) 



C 



g/CSs(^' ^' /^i) 



2n 



(1 



x] 



6{x 



7T 

4 

2 



(33) 



Comparing with the results in |7|, we find that we can reproduce the A^l^, A^SS' -^css' 
Cpgg, which are all the coefficients and functions needed for resummation at NLL order. 

(2) 

-^css ^^^^ ^® needed to resum NNLL. However, to fully achieve NNLL resummation we need 



.(2) 



?(1) 
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to calculate the matching coefficient at the lower scale up to order in SCET [l9|. We 
leave this to a future publication. 

Following the above, the resummation for SIDIS can be performed similarly. As we 
stated earlier, the anomalous dimension will be the same. The only difference is the process- 
dependent matching coefficients at and /i| in the resummation formula. For DIS, one 
has the one-loop result. 



Zn 



a.9 



4-^ 
12 



(34) 



which leads to the result for the C 



(1) 

g/CSS 



2ti 



in DIS, 

[(1 - x) - 4(5(x - 1)] . 



(35) 



These results agree with those from the conventional resummation approach [20| • 



III. STANDARD-MODEL HIGGS PRODUCTION 

Transverse-momentum dependence of the Standard Model Higgs production can also be 
studied through resummation of large double logarithms. Higgs production, for a large range 
of Higgs mass, can be described by an effective action with a pointlike coupling between the 
Higgs particle and gluon fields. The effective coupling is of course scale dependent, bal- 
ancing the renormalization dependence of the composite operator. In general, this effective 
lagrangian can be written as |lq . 

Cnaa = CEw{Mt)CT{Mt, /i)i/F;,F'^^^ , (36) 

where H is the scalar Higgs field and F^j, is the gluon field strength tensor. Cew represents 
the electroweak coupling coefficient from the heavy-top-quark loop calculation, while Ct 
comes from the strong interaction radiative corrections. The coefficient Ct will depend, in 
general, on the top quark mass and the renormalization scale /i. To our interest, we quote 
at two-loop |2l[22j. 



4tt [ An 



(47r)^ 



In ^ (7Cl - IIC^C^ + mfTpCF) + 



(37) 



Here we have omitted the constant terms of order because they do not contribute to the 
renormalization group running at two-loop order. 

In SCET, the Higgs production cross section can be calculated from its coupling with the 
collinear gluon fields, 

Ch-scet = C{Mh, /i)// J-„V^"*^^" , (38) 

where T'^^^y represent the n and n collinear gluon field strength tensors in SCET 
C{MH,fi) is the matching coefficient which contains the coupling of Higgs boson to gluons 
in full QCD, 

C{MH,fi) = CEw{Mt)CT{Mt,fi)CG{MH,fi) . (39) 
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The last factor comes from matching between operator HF^j^F"'^'^ in full QCD and 
H^nfj.u^"'"'^'^ SCET. Because Cew has no QCD effects, we will not discuss it further 
in this paper. Ct and Cg contain the QCD evolution effects, and thus the anomalous 
dimension of the SCET operator HJ^^^^J^"-'^^'^ is the sum of the two: 

7i = 7t + 7g , (40) 

where jt and •jg are defined as 



It 



d In fi 



7G = —7] • (41 



From Eq. (jSTj), it is easy to show that 



7r = -2/3o-^ + j-^ [-213, - 2/3o(5Ca - 3^,.) 

+2{7Cl-llCACF + 8NfTFCF)] , (42) 

up to two-loop order. To calculate 7g, we follow the calculation for the Drell-Yan process in 
the previous section. Using the gluon form factor in 0,11^, we get the relevant anomalous 
dimension, 



7g = — CaIh^ 



+ (4^ 



(43) 



where fgl'' can be obtained from the above fj^l^ by color- factor exchange {Ca ^ Cp) jl6l |. 
The anomalous dimension 7^3 



^(1) 
'99 



Cl (^y + 24C(3)^ - y iV/TpC^ - SNfTFCp , (44) 



is proportional to 5{1 — x) in the gluon splitting function. (3i is the two-loop beta function 
defined before. 

Including the matching at lower scale /i^, the result for Higgs production W{h, Mfj, xi, X2) 
at low transverse momentum can be written as, 

W{Mjj, h) = CUmI, a,{Ml))e^^ {C, ® g) (xi, 6, /ii) (C, ® g) (x2, 6, A) , (45) 

where Cew and the leading factor in Ct have been absorbed in the Born cross section, and 
the remainder Ch is 

Ch{MI, as{Ml)) = 1 + [Wa - ^Cf + Ca (C(2) + vr^)] , (46) 
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up to one-loop order. The Sudakov suppression form factor has the same form as that for 
DY, and the expansions of the A and B are, 





= Ca 






Am 


= CaK 




1 - 
~ ~16 - 



7^^ + + 2/^0 (5Ca - 3Cp) - 2(3, 
-2 {7Cl - 1 \CaCf + 'iN^TpCp) ] (47) 

With these results we can reproduce the conventional resummation for Higgs-boson produc- 
tion at NLL order [l9| . 

The coefficient Cg to one-loop order is calculated from the matching at the scale of /i^, 
with real and virtual contributions. It is needed for resummation at NLL. The result is 

X, ^l) = 6{1 -x) + ^6{x - 1) (^-^) , (48) 

where we have chosen hl = Ci/b with Ci = 2e^'"^ to eliminate the large logarithms. The 
corresponding coefficient for CSS resummation is 

C,/css(&, X, /xi) = 6{1 -x) + ^5(1 - x) [(5 + n')CA - 3Cp] (49) 
which is consistent with |l 



IV. CONCLUSION 

In this paper, we demonstrated how to perform resummation of large double logarithms 
in hard processes involving low transverse momentum in the framework of effective field 
theory. The newly-discovered theory, SCET, is naturally suited for various hard processes 
with two or more well-separated momentum scales . As an example, we studied the Drell- 
Yan resummation in detail, and outlined how to extend the procedure to other processes as 
well. We calculated the relevant anomalous dimensions of the effective current up to two- 
loop order, and reproduced the conventional resummation results to NLL accuracy. This 
method can be certainly extended to higher order as well, for example, NNLL. However, 
this requires calculation of the matching coefficient at the intermediate scale within SCET 
up to two-loop order. We leave these extensions for a future publication. 
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